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' ^; Abstract 
^ ■ 

^ ' We consider a generalized discrete cj) model and demonstrate that it can support exact moving 

^ : 

. kink solutions in the form of tanh with an arbitrarily large velocity. The constructed exact moving 

m ■ 

\ solutions are dependent on the specific value of the propagation velocity. We demonstrate that 

(N ■ 

' in this class of models, given a specific velocity, the problem of finding the exact moving solution 

oo : 

■ is integrable. Namely, this problem originally expressed as a three-point map can be reduced to 
• ^ , a two-point map, from which the exact moving solutions can be derived iteratively. It was also 

X ■ 

H ' found that these high-speed kinks can be stable and robust against perturbations introduced in 

the initial conditions. 
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I. INTRODUCTION 



Solitary waves play important role in a great variety of applications because they are 
robust against perturbations and they can transport various physical quantities such as mass, 
energy, momentum, electrical charge, and also information l|, l2|. Many popular continuum 
nonlinear equations support traveling solutions, but for their discrete analogs the existence 
of traveling solutions was not systematically studied until recently. Previously, only a few 
integrable lattices were known to support moving solitons, e.g., the Toda lattice [3] and 
the Ablowitz-Ladik lattice [4], the former being the discrete version of the KdV equation 
and the latter of the nonlinear Schrodinger (NLS) equation. While integrable lattices of 
other key continuum models are also available [5|, |6|], these integrable situations are still 
rather exceptional and typically not directly relevant to experimental settings. On the other 
hand, recently several wide classes of discrete models supporting translationally invariant 
(TI) static or stationary solutions have been discovered and investigated for the discrete 



Klein-Gordon equation 
NLS equation |2o|, bj 
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Static solitary waves in such lattices possess the 
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22J, which means that the solitary waves 



moving with vanishing velocity can be regarded as the exact solutions to these discrete 



equations (although there are issues for finite but small velocities as explained e.g., in 



27l|). Moving solutions propagating along a lattice at finite velocity have also been found 
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and with the help of specially tuned numerical approaches in 24, 



32|. Traveling bright solitons in the NLS lattice model with saturable 



nonlinearity are also very interesting 



21, 
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34| . For this model the existence of soliton 



frequencies with vanishing Peierls-Nabarro energy barrier was demonstrated. It was shown 
that for the specific frequencies the static version of the model equation coincides with 



that of the integrable Ablowitz-Ladik equation [3J]. This implies in the mapping analysis 
a possibility to reduce the three-point map to a two-point map. However, the studied 
model equation is not integrable and while the solitary waves moving with small velocity 
were numerically obtained for the specific frequencies, the question about the existence of 
solitary waves moving with finite velocities remains open. 

Existence of lattices supporting the static Tl solutions and the exact solutions moving 
with finite velocity poses a natural question, about whether there is a relation between them. 
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In the present study we answer this question in the positive, creating a hnk between these 
two classes of solutions below. 

On the other hand, as regards the solitary wave motion, another interesting question is 
the existence of upper bound for the propagation velocity. In view of Lorentz invariance, 
clearly the solitary waves in continuum nonlinear equations of the Klein-Gordon type have 
limitations on their propagation velocity. One of the questions that we examine below is 
whether a similar restriction exists in the corresponding discrete models. Since the Lorentz 
invariance is no longer a symmetry of the discrete nonlinear systems, strictly speaking there 
need not be any restriction on the propagation velocity. However, naively one would expect 
that typically the propagation velocity in the discrete nonlinear systems would be similar 
to the ones supported by their continuous counterparts and not develop supersonic values. 
Recently, three of us [22l] confirmed this naive expectation in a generalized discrete NLS 
equation. A natural question, however, is whether this always holds true. Although in some 



case examples, this has already been addressed [26|], here we address this question (and its 
negative answer) more systematically in the framework of the recently studied discrete (p^ 
model Q: 



A: 



-^0n+10n-l((/'n+l + (pn-l) " ^(0^+1 + ^Ll). (l) 

with the model parameters satisfying the normalization constraint 

6 

J2Ak = X. (2) 

k=l 

In Eq. ([1]), (f)n{t) is the unknown function defined on the lattice x„ = hn with the lattice 
spacing h > 0. Without loss of generality it is sufficient to consider the cases A = 1 and 
A = -1. 

Our results will be displayed as follows. Analytical results are collected in Sec. [Tll 
Examples of the exact moving solutions to Eq. ([T]) are given in Sec. Ill Al to Sec. Ill CI 
Then, in Sec. IIIDt we derive a continuum analog of Eq. ([1]) and find its kink solution. In 
Sec. Ill El spectrum of the vacuum solution of Eq. ([1]) is obtained. The relation between 



the exact moving solutions and integrable maps is established in Sec. IIIFI Section |TTT] is 
devoted to the numerical analysis of stability and robustness of moving kinks. We present 
our conclusions in Sec. [IVl 



II. ANALYTICAL RESULTS 

Our approach to deriving analytical solutions will be, at least in the beginning of this 
section somewhat similar in spirit to the inverse method of [sS^. That is, we will postulate 
the desired solution and will identify the model parameters for which this solution will be a 
valid one, as we explain in more detail below. 

A. Moving sn solution 

It is easy to show that an exact solution to Eq. ([T]) is 

(pn = Asia[P{hn + hxo — vt),m] , (3) 

provided the following relations are satisfied: 

A5ias{2hi3, m) = —A^ias^hp, m) , 

A^A^ = -2m(3\^ , A = 0, (4) 

2 o 
— — = 2Ains(hl3, m)ns(2/i/5, m) + A2ns^(/i/5, m) 
A'^h^ 

—Ascs^hp, m)ds(/i/5, m) 
-A5[cs(2/i/3, m)ds{2hp, m) - ns'^{2hp, m)] , (5) 

2 — A/i^ 

— = A4ns^(hp, m) + yl2Cs(/i/9, m)ds(/i/5, m) 
A-^h^ 

-A3ns'^{h/3,m) + {l + m)-^. (6) 

In Eq. ([2]), A is the amplitude, v is the propagation velocity, /? is a parameter related 
to the (inverse) width, xq is an arbitrary position shift, < m < 1 is the Jacobi elliptic 
function (JEF) modulus. Further, ns(x, m) = l/sn(x, m), cs(x, m) = cn(x, m)/sn(x, m) 
and ds{x,m) = dn(x, m)/sn(x, m), where sn(x, m), cn(x, m), and dn(x, m) denote standard 
JEFs. 
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Solution parameters xq and m can be chosen arbitrarily. Equations (jl]) to ([6]) establish 
five constraints from which one can find the three solution parameters A, v, and (3, and two 
constraints on model parameters h, X, and Ai {i = 1, 2, 6). 

It is possible to construct some other moving JEF solutions, for example, moving cn and 
dn solutions (and hence hyperbolic pulse solutions of the sech type) but we do not discuss 
these here. 



B. Exact moving kink solution 

In the limit m — 1, the above moving sn solution reduces to the moving kink solution 

0(x, t) = tanh [P {hn + hx^ — vt)] , (7) 

and the relations (jl]) to ([6]) take the simpler form 

/?2 = ^, 2A, = -A,{l + T), A = 0, (8) 
2T 

— = A4(l + r) + A2 + ^(l + 2T-T2), (9) 

^l:iJ^^ = A, + A,{1 -T) + ^(1 + T) - A,T , (10) 



where 



T = tanh2(/i/5). (11) 



One can see that the solution is defined onl y i 



Ai < 0. 



It is worth pointing out that the static 22| and the moving JEF as well as hyperbolic 



soliton solutions exist in this model in the following seven cases: (i) only A2 nonzero (ii) only 
nonzero (iii) ^3, Ar, nonzero (iv) A2, A4 nonzero (v) A2, A3, A^ nonzero (vi) ^3, A4, A5 
nonzero (vii) as discussed above A2, A^, At^^ A^ all nonzero. Further, while in the static case, 
solutions exist only if Ai, = 0, in the moving case solutions exist only if A^ = 0, 7^ 0. 
It may be noted that these conclusions are valid for sn, cn, dn as well as tanh, sech solutions. 
Moreover, while in the static case, the kink solution exists only if A > 0, the moving kink 
solution is possible even when A < 0. In fact, it turns out that in cases (i), (ii) and (iii) kink 
solution with a large velocity v is possible only if A < 0. 

While Eqs. (IE]) to (fTOl) give a general set of restrictions on the model parameters sup- 
porting the exact moving kink, below we extract and analyze three particular cases where 
the restrictions attain a very simple and transparent form. 



Case I. Only Ai, A2, and A4 are nonzero. These parameters are subject to the constraint 
of Eq. ([2]) and one can take Ai < and A4 < Ai ~ \ + 2/h^ as free parameters. From 
equations ([BD to (fTU]) and we get 

P' = ^^ A, = ^-^, A, = X-M-M. (12) 

We note that the kink velocity v can have an arbitrary value in case A > 0. 

Indeed, for chosen A and /i, one can take any v and then, using the expressions of Eq. 
f|T2|) . find subsequently the inverse kink width j3 and the remaining model parameters 
and A2. 

As a subcase of case I, one can take only Ai and A2 nonzero. From Eq. ([T^ we get 

/32 = ^, 2T = h^\ + 2v^ei3\ A2 = \-A^. (13) 

From here it is clear that the kink solution with large v is only possible in this case if A < 0. 
Similar arguments are also valid in cases (ii) and (iii) discussed above. 
Case II. Another interesting case is when {A2, A^, nonzero): 

2 2 

^2 = A- — -Ai, A3 = X- — - Ai + Ai, 

A, = -\ + ^ + A^- 2A4, Ae = 0, (14) 

with Ai and A4 being free parameters. It is clear that the constraint of Eq. ([2]) is satisfied 
for any Ai and A4. Conditions ([8]) to (fTOj) reduce to 

-Ai _A,-\ 

Thus, we have a two-parameter set of moving kinks. Kink solution exists for Ai < and 
A4 > 2/h?. In this case too the kink velocity v can obtain any value. 



= A = (15) 



As it was found in 19|, for 

= 471A, A2 = 672A, A3 = (1 - 471 - 872)A, 

A4 = A5 = 0, Ae = 272A , (16) 

with arbitrary 71 and 72, the model Eq. ([T]) has the Hamiltonian 

rr -S^l'^n . (0n - 0«-l)^ A A 3 , ^ , 4 
n 

+72</'n</'n-l(0n + 0n-l) +(^^71 - 272 j ^^^Ll , (17) 



and hence energy is conserved in this model. Simple analysis of Eqs. ([8]) to (fTOl) suggests 
that among the models supporting the exact moving kink solutions there are no Hamiltonian 
models. 



C. Exact moving trigonometric solution 
and moving four-periodic solution 

Exact moving trigonometric solution. The discrete model of Eq. ([T]) supports an exact 
moving solution of the form 

(j)n = A sm[hp{n + Xq) — vt] , (18) 
even when Ai, i = 1, 6 are all nonzero provided 

-2 + 2 cos{hp) + 

= A^h^ sm\h(3)[A3 - ^4 + {3Aq - A^) cos{h(3)] , (19) 

Ai + Ai + A3 cos(2/i/?) + {A2 + A5) cos{hp) 

+Aq cos{hp) [4 cos^hp) - 3] = . (20) 

From here, one can work out specific relations in the case of various models. For example, 
consider the Hamiltonian model with the parameters satisfying Eq. flTBl) . It is easy to check 
that the moving sine solution exists in this model provided 

-2 + 2cos{hp) + 

= A^Xh^ sin^ihp) [1 - 471 - 872 + 672 cos(/i/?)] , (21) 

471 + (1 - 471 - 872) cos{2hf3) + 872 cos3(/i/?) = 0. (22) 

Exact moving four-site periodic solution. For jSh = 7t/2 the above moving trigonometric 
solution reduces to the moving four-periodic solution of the form 

'IT 



')n = A sin 



-{n + xq) - vt 



(23) 



with 

Aih^ 



+ = A3 , A^ = — . (24) 



For the Hamiltonian model, Eq. (12^ reduces to 

71 + 72 = ^ , A^= i;^^ . (25) 



1 A/i2-2 + t;2 

47iA/i2 



Even in the Hamiltonian case the trigonometric solution can have an arbitrarily large veloc- 
ity. 

D. Approximate moving kink solution 

The discrete model of Eq. (jl]) with finite h can support the solutions varying slowly with 
n (the same class of solutions can be studied near the continuum limit, using /i as a small 
expansion parameter). For such solutions one can use the expansion 0(/i) cj) + hcf)^ + 
{l/2)h'^(l)xx and substitute in Eq. ([T]) 0„ = 0(0), = (p{±h) to obtain 

<Ptt = <Pxx + A0 - X(l)'' - hi''B<p''<P,, - h^D(t)(t)l , (26) 

where 

B = A2 + 2A3 + 2A4 + 3A5 + 3Ae , 

D = A3-A,-A, + 3Ae, (27) 

and Eq. ([2]) was used. For B = D = Eq. fl26l) reduces to the continuum (j)^ equation. This 
equation also stems from Eq. (1261) in the limit h ^ 0. 
Equation (l26l) supports the following kink solution 

0(x,t) =tanh[/?(x + Xo-ft)], (28) 
provided the following two conditions are satisfied 

^'- 2(l-.') + ft^C - ^ + " = 

This is the quasi-continuum analog of the solution of Eq. ([7]). Let us analyze the case 
of A = 1, which corresponds to the double-well potential. For D = B = or/and h = 
we have the classical continuum 0'^ kink with the propagation velocity limited as |f | < 1. 
On the other hand, for nonzero D = —B and nonzero h, there is no limitation on the kink 
propagation velocity. In other words, fast kinks are possible in Eq. (!26l) only in the presence 
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of several competing nonlinear terms. Note that in case B,D ^ 0, Eq. (l26l) is not Lorentz 
invariant. 

For wide kinks, i.e. when (3 ^ h, one can use the kink solution of Eq. fl26l) to write 
down the approximate solution to the discrete model Eq. ([T]) of the form of Eq. ([7]), whose 
parameters are given by Eq. fl29p . The corresponding wide, fast kinks are possible in Eq. 
([1]) only for finite h and only ioi D = —B ^ Q. 

Substituting the Hamiltonian conditions of Eq. f|T6|) into Eq. fl27j) we find B = 2D while, 
according to Eq. (!29l) . the approximate moving kink solution exists at -B = —D, i.e., it 
exists only in non-Hamiltonian discrete models. 

E. Spectrum of vacuum 

The discrete model of Eq. ([T]) supports the vacuum solutions 0„ = ±1. The spectrum 
of the small- amplitude waves (phonons) of the form ~ exp{ikn ± iut), propagating in 
the vacuum, is 



where k denotes wavenumber, u is frequency and B is as given by Eq. fl27|) . 

Another vacuum solution, 0„ = 0, supports the phonons with the dispersion relation 



and this vacuum is stable for A < 0. 
F. Relation to TI models 

Note that for t> = the exact moving solutions Eq. (j3]), Eq. ([71), and Eq. ( |T8|) reduce to 
the static TI solutions, i.e., solutions which, due to the presence of arbitrary shift xq, can 
be placed anywhere with respect to the lattice. TI solutions of Eq. ([T]) were discussed in 
detail in [l^. Particularly, it was established that any TI static solution can be obtained 
iteratively from a two-point nonlinear map. Let us demonstrate that the exact moving 
solutions (including JEF solutions) can also be derived from a map. 

Let us generalize the ansatz Eq. ([3]) and look for moving solutions to Eq. ([T]) of the form 




(30) 




(31) 



= AF{^) , e = P{hn + hxo - vt) , 



(32) 
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with constant amplitude A, where F is such a function that 

(f(t>r, 



-(3'v\ao(j)n + ai(Pi), (33) 



with constant coefficients ao, «!• If this is the case, the substitution of Eq. (!32|) into Eq. 
([1]) will result in a static problem essentially identical to the static form of Eq. ([1]), namely, 
in the problem 

2 

= J^i^n-l - 20„ + (pn+l) + 2(A + aoP^V^)^n 
-2(Ai - aiP^^m - A24>l{<Pn-l + 0n+l) 
-A'i(t)n{(t)n-1 + 4>n+l) ~ 2A40„_i0„0„+i 
-A50„„10„+l(0„_l + - ^6(0^.1 + 0^+1 ). (34) 

This effective separation of the spatial and temporal part with each of them satisfying 



appropriate conditions is reminiscent of the method of 36 1. 

The Jacobi elliptic functions and some of their complexes do possess this property. For 
example, Eq. (1331) is valid for 

2m 

tto = 1 + m, tti = - — 

2m 

ao = 1 - 2m , tti = —TT 

A"^ 

2 

ao = m - 2 , ai = 

2 

ao = 1 + m, ai = - — 
2fm- 1) 

ao = 1 — 2m , ai = 



for 


F = 


sn; 


for 


F = 


cn; 


for 


F = 


dn; 


for 


F = 


1 _ 

sn' 


for 
for 


F = 
F = 


1 

cn' 
sndn 



A^ ' 
2 

ao = 2(2m - 1) , ai = — - , for F = (35) 

A^ cn 

In the limit m = 1 we get from the first line of Eq. fl35l) 

ao = 2 , ai = —2 , for F = tanh . (36) 

A solution moving along the chain continuously passes through all configurations between 
the on-site and the inter-site ones. If one finds a static solution to Eq. (IMl) that exists for 
any location with respect to the lattice, then one finds the corresponding moving solution to 
Eq. ([1]). Various static solutions to Eq. flM|) that have the desired property of translational 
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invariance have been recently constructed. Let us discuss the two simple cases specified in 
Sec. HTBl 
Case I. 

For nonzero Ai, A2, and A4 = X — Ai — A2, Eq. ([T]) reduces to 
d^(h„, 1 



- 2(j)n + 0„+l) + X(f)n " Ai( 



while from Eq. (IMl) we get 

1 

h 



Setting 



we reduce Eq. 



with 



2 y^n-l - '2(j)n + (j)n+l) + (A + Oof^'^v'^) 



-{Ai - ai(3^v 



2„,2n i3 



0n(0n+l + 0n-l) - A0n-10n0n.+l = . 



2„.2 



Ai = aiP V 



ao = — cti 



to the integrable static equation 



17| 



-(P- A 



(37) 



(38) 



(39) 



(40) 



P = \-Ai. 



(41) 



Any static TI solution to Eq. ( HOj) generates the corresponding moving solution to Eq. (1371) . 
provided that Eq. fl33l) is satisfied. 



All static solutions of Eq. (140|) can be found from its first integral 19| 



^2 .2 _ y^n^Pl+l 



CY 



(2 - P/i2)2 - C/i4(p - A2f 



2{2 - Ph^) + Ch^A2{P - A2) 
Y = h\P-A2 + A2Z), 



(42) 



where C is the integration constant. The first integral can be viewed as a nonlinear map 
from which a particular solution can be found iteratively for an admissible initial value (po 
and for a chosen value of C. 
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Thus, Eq. (H2I) is a general solution to Eq. (HOj) , while to make it also a solution of Eq. (138|) 
one needs to satisfy Eq. fl39l) . This can be achieved by calculating (^(pn/d^'^ from the map 
Eq. and equalizing it to — ao0n+i "C^i^n+i according to Eq. (15^ . To calculate (P(j)n/d^^ 
we consider a static solution centered at the lattice point with the number n so that = 0. 
From the map Eq. fH21) one can find 0„+i(O) = a/ CF/ (2 — P/i^). We then consider the initial 
values (pn = —£ and 0„ = e and similarly, from the map Eq. (H2l) . find the corresponding 
values of 0„+i(— £:) and 0„+i(£:). Finally, we calculate the second derivative at the point 
n + 1 as d^(t)n/d(^ = lim (-e) - 20„+i (0) + 0„+i (e:)] /e:^ = -ao0n+i - where 
the values of the coefficients are 

CY' - (2 - Ph'f {Z^ - 1) 2 

= (2 - p/.^) ' "^ = -c- ^^^^ 

Let us summarize our findings. From Eq. fH5]) and the second expression of Eq. (1391) . 
taking into account Eq. fHTl) . we find the relation between the integration constant C and 
model parameters Ai, y42. A, and /i when a moving solution is possible: 

{I - C) P {2 - PK'Y 
X [/^^(C - l)p3 + ^/i4^^ ^5^^p _ _ 2^2^ 

+4 (2/I2 - 3P - AaP/i^ + 2p2/i2) 1=0. (44) 

The solution profile is found for this C from Eq. fH2]) iteratively for an admissible initial 
value 00, and the propagation velocity v is found from the first expression of Eq. fl39l) and 
the second expression of Eq. (H3l) . Following this way, any moving solution to Eq. (!37|l 
can be constructed (possibly, except for some very special solutions that may arise from 



factorized equations [19|]). Note that in [17| we could express many but not all the solutions 
of Eq. fHU]) in terms of JEF. The approach developed in this section allows one to obtain 
iteratively even those moving solutions whose corresponding static problems were not solved 
in terms of JEF. 

The simplest case is C = 1 when Eq. fj44|) is satisfied and this case corresponds to the 
kink solution. From Eq. (H3!) we find 

ao = 2, «! = -2 , (45) 

which coincides with Eq. (1361) . The map Eq. fH2|) in this case reduces to 

± - 0n+l) - 0n0n+l + 1 = 0. (46) 
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One can check that Eq. (H6|) supports the static solution 0„ = tanh[/i/3(n + xq)] with 



and this coincides with the second expression of Eq. f[T2|) . The static solution Eq. (j4j 
that can also be found iteratively from Eq. fl46l) . gives the profile of the moving kink that 
satisfies Eq. f l57|) . The kink velocity v is found from Eq. fl5^ and Eq. fH5]) and this agrees 
with the first expression of Eq. f|T2l) . 
Case //. 

The model parameters are related by Eq. f|T4|) with and being free parameters. 

For Ai = -2(3^v^, Eq. §^ reduces to 

2 

^ (0„_i - 20„ + 0„+i) + 2 (A - Ai) (pn 

-A24>1 (0n-l + 0n+l) - ^30n (0^1 + 0'+l) 

-2A40„_i0„0„+i - A50„_i0„+i(0„_i + 0„+i) = 0, (48) 

I 1 

which is a particular form of the case (vii) static equation studied in [19]. The first integral 
of the static model (vii) for the general case is not known, but it is known for any particular 
JEF solution [l9|. Let us further simplify the problem considering only the kink solution of 
Eq. (I48l) . for which one has 

tanh2(/i/5) = . (49) 

This relation coincides with the second expression of Eq. ( ITSi) . From the known static tanh 
solution one can deduce the following two-point map that generates this solution for any 
initial value |0o| < 1- 

0„ ± v/(Ai - A)M5 



■>n+l 



(50) 



where one can interchange 0„ and 0„+i and take either the upper or the lower sign. 

The static kink solution to Eq. ( l48l) with jS satisfying Eq. (H9l) . that can also be found 
iteratively from Eq. (!50|) . gives the profile of the moving kink that satisfies Eq. ([1]) with the 
parameters related by Eq. f|T^ . The kink velocity v is found from Eq. fl39l) and Eq. fHSl) . 



III. NUMERICS 



Let us analyze the kink solutions for the cases I and II described in Sec. Ill B[ Only the 
case of A = 1 will be analyzed. 
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FIG. 1: Isolines of equal kink velocity on the plane of model parameters Ai, for the Case I at 
\ = h = l, A2 = \- Ai- Ai, A-i = A^ = Aq = {). 




-15 -10 -5 5 10 15 



1.7 1 ■ ' ■ ' ' ■ 1 

0.0 0.5 1.0 1.5 2.0 

t 



FIG. 2: Kink velocity as a function of time for the Case I. The exact kink velocity for chosen model 
parameters is f = 10 but for setting the initial conditions we put in Eq. ^ v = 9.8 (dots) and 
V = 10.2 (open circles). In both cases, regardless of the sign of perturbation, the kink velocity 
rapidly approaches the exact value. The inset shows the moving kink profile. Model parameters: 
\ = h = l, Ai = -20, Aa ^ -222.1, A2 = \- Ai- A^, A3 = A5 = Aq = 0. 
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In our simulations we solve the set of equations of motion, Eq. ([T]), numerically with a 
sufficiently small time step r using the Stormer integration scheme of order O(r^). Initial 
conditions are set by utilizing Eq. ([7]) with various parameters. Anti-periodic or fixed 
boundary conditions are employed. 

A. Exact moving kinks in Case I 

In this case, as it was already mentioned, the exact moving kink solution exists for the 
free model parameters satisfying Ai < and ^4 < Ai — X + 2/h'^ (see Fig. [T] where for 
X = h = 1 we show the isolines of equal kink velocity on the plane of model parameters Ai, 
A4). On the line Ai = 0, according to the first expression of Eq. (fT2|) . the kink velocity v 
vanishes. On the line A4 = Ai — X + 2/h'^ the kink velocity also vanishes. This is so because, 
as it can be seen from the second expression of Eq. fll2l) . on this line we have T = 1, i.e. 
/3 — s> 00 and kink width vanishes. On the line A4 = Ai — X + we also have B = 
which corresponds to vanishing of the width of the phonon spectrum given by Eq. ( |30ll . 

The vacuum 0„ = ±1 is stable, i.e. the spectrum Eq. (!30|) does not have imaginary 
frequencies, if A4 < Ai + 2/h?, i.e., it is stable in the whole region where the exact moving 
kink solution is defined. 

As it was already mentioned, the kink propagation velocity is unlimited and from Fig. 
[Done can see that 1^441 increases for higher kink velocities while Ai can have any negative 
value. For example, for f = 10 the moving kink exists for Ai <Q and A^ < —221.9. 

Let us now turn to the discussion of the stability of moving kinks. 

Case of X = h = 1 and v = 10. 

Moving kink solutions are stable not in the whole range of parameters of their existence 
but for each studied value of propagation velocity we were able to find a range of parameters 
where the moving kink is stable in the sense that even a perturbed kink solution (with a 
reasonably large perturbation amplitude) in course of time tends to the exact solution (i.e., 
it is effectively "attractive"). The evolution of the kink velocity in such self- regulated kink 
dynamics is shown in Fig. [2] for a velocity as large as f = 10 for Ai = —20 at X = h = 1. 
For this choice, we have from Eq. ([12]) /3 ^ 0.316, ^4 ^ -222.1, A2 ^ 243.1. The moving 
kink profile is shown in the inset of Fig. [2J The perturbation was introduced into the 
initial conditions by setting in the exact solution of Eq. (JTj), a "wrong" propagation velocity. 
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smaller or higher than the exact value v = 10. It can be seen that in course of time the 
propagation velocity approaches the exact value regardless of the sign of perturbation. The 
increase of kink velocity launched with v = 9.8 may look counterintuitive but one should 
keep in mind that moving kinks are the solutions to a non-Hamiltonian (open) system with 
the possibility to have energy exchange with the surroundings with gain or loss, depending 
on the trajectories of particles (see, e.g., |10|). 

We found that the absolutely stable, self-regulated motion of kink, similar to that shown 
in Fig. [2], for v = 10 and X = h = 1 takes place within the range of —164 < Ai < —1.1. The 
inverse kink width at the lower edge of the stability window is /? ~ 0.906, which corresponds 
to a rather sharp kink, while for Ai = —1.1 (the upper edge of the stability window) one 
has (3 ~ 0.0742 and the kink is much wider than the lattice spacing h. For comparison, the 
kink shown in the inset of Fig. [2] has (3 ~ 0.316. For Ai < —164 the moving kink solution 
becomes unstable and displacements of particles behind the kink grow rapidly with time 
resulting in the stopping of numerical run due to floating point overflow. For — 1.1 < Ai < 
kink dynamics is as described in Sec. 6 of 1^ for the non-Hamiltonian case. In this region 
of parameter Ai, after a transition period, the kink starts to excite in the vacuum in its wake 
a wave with constant amplitude. In this regime, the kink attains a constant velocity whose 
value is, generally speaking, different from that prescribed by Eq. (fT2!) . Kink dynamics in 
the two unstable regimes described above will be illustrated below for the kinks moving with 
V < 1. 

Case of X = h = 1 and v < 1. 

Similar results were obtained for the kinks moving with small velocities {v < 1). For 
smaller velocities the range of Ai with stable, self-regulated motion becomes narrower. For 
instance, a kink with v = 0.8 is stable (in the above mentioned sense) for —2.4 < Ai < —0.10, 
while the one with v = 0.5 for -0.19 < Ai < -0.06. 

In Fig. [3] we show the time variation of the moving kink profile to demonstrate the 
instability of the exact kink solution moving with v = 0.8 at Ai = —2.6. This solution is 
unstable and, due to the presence of rounding error perturbation, the wave behind the kink 
is excited. The amplitude of the wave rapidly grows with time and it becomes noticeable in 
the scale of the figure at t > 100. At t ~ 130 the numerical run stops due to the floating 
point overflow. Kink velocity is nearly equal to f = 0.8 practically until the collapse of the 
wave behind it. 
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In Fig. m we show the kink velocity as a function of time to illustrate the instability of 
the exact moving kink solution at f = 0.8 and Ai = —0.07. This solution is unstable and, 
due to the presence of rounding error perturbation, at t ~ 400 it starts to transform into 
an oscillating kink (with the period T ^ 5.01) moving with different velocity and exciting a 
constant-amplitude wave behind it. The transformation is essentially complete by t ~ 5000 
and the kink velocity becomes v ~ 0.96. The inset shows the profiles of an oscillating 
moving kink in the two configurations with the most deviation from the average in time 
configuration. 

Case of X = 1, h = 0.1, and v = 10. 

Stable, self-regulated motion of high-speed kinks was also observed for as small lattice 
spacing as h = 0.1. This may appear surprising at first sight because for small h one would 
expect the discrete model to be close to the continuum 0^ model where propagation with 
the velocity faster than = 1 is impossible. But looking at Eq. (126!) . one can notice that 
the two last non-Lorentz-invariant terms have the coefficients h^B and h^D, and they are 
not small even for small h if B and D are large. The high-speed kink in the considered 
case indeed exists for B ^ —19833.1 and D ^ 20034.1. These values correspond to the 
following parameters considered in this numerical run: Ai = —200.0, ^4 ^ —20034.1, 
A2 = X-Ai-A^^ 20235.1, A3 = A5 = Aq = 0, p = 1. 

B. Exact moving kinks in Case II 

In the case II, the model parameters are related by Eq. f|T^ and we find B = —2A4^+Q/h?. 
As it was already mentioned, the exact moving kink solution exists for Ai < and A^ > 2/h'^. 
Isolines of equal kink velocity on the plane of the free model parameters Ai, A4 can be seen 
in Fig. [5]for X = h = 1. In the Case II, similarly to the Case I, on the borders of the existence 
of the kink solution the kink velocity vanishes. Moreover, on the border A4 = 2/h'^, both 
the width of the phonon spectrum given by Eq. fl30l) and the kink width vanish, and this is 
also similar to what we saw in Case I. 

As was mentioned above, the exact moving kink solution in the Case II can also have 
arbitrary velocity. Kink dynamics in the Case II was observed to be qualitatively similar to 
that of the Case I. For the kink with speed as large as f = 10 for X = h = Iwe found that the 
absolutely stable self-regulated motion is observed within the range of —82 < Ai < —1.0. 
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FIG. 3: Change of the moving kink profile demonstrating the instabihty of the exact solution for the 
following set of model parameters: X = h = 1, Ai = —2.6, A4 ~ —2.538, A2 = X — Ai — A4 ~ 6.138, 
A^ = A^ = Aq = 0. Exact kink parameters are v = 0.8, /? ~ 1.425. This solution is unstable 
and, due to the presence of rounding error perturbation, the wave behind the kink is excited. The 
amplitude of the wave rapidly growths with time and it becomes noticeable in the scale of the 
figure at t > 100. At t ~ 130 the numerical run stops due to the floating point overflow. 

The corresponding inverse kink width varies in the range 0.640 < (3 < 0.0707. For Ai < —82 
the moving kink solution becomes unstable and displacements of particles in the wake of the 
kink grow rapidly with time resulting in the termination of numerical run due to floating 
point overflow. For —1.0 < Ai < 0, after a transition period, the kink starts to excite in the 
vacuum behind itself a wave with constant amplitude. In this regime, the kink propagates 
with a constant velocity whose value is, generally speaking, different from that prescribed 
byEq. (USD. 

C. Exact moving four-periodic solution 

We have studied the dynamics of the exact moving four-site periodic solution Eq. (l23ll 
in the Hamiltonian model, i.e., with the parameters satisfying Eq. (125|) . Periodic boundary 
conditions were employed for a chain of 40 particles. A small perturbation was introduced 
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FIG. 4: Kink velocity as a function of time demonstrating the instability of the exact moving 
kink solution for the following set of model parameters: \ = h = 1, Ai = —0.07, ^4 ~ —18.28, 
A2 = \ - Ai- A^K. 19.35, A^ = A^ = A^ = 0. Exact kink parameters are v = 0.8, j3 0.234. 
This solution is unstable and, due to the presence of rounding error perturbation, at t ~ 400 it 
starts to transform to an oscillating moving kink (with the period T ~ 5.01). The transformation 
is essentially complete by t ~ 5000. Inset shows the profiles of the oscillating moving kink in the 
two configurations with the most deviation from the average in time configuration. 




FIG. 5: Isolines of equal kink velocity on the plane of model parameters Ai, ^4 at A = /i = 1 for 
the case 11 when model parameters satisfy Eq. (I14p . 
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in the amplitude of particles at t = and we always observed a growth of deviation of the 
perturbed solution away from the exact one. We varied the free model parameter 71 and 
the solution parameter v in wide ranges at A = /i = 1 but could not find a stable regime. 



IV. CONCLUSIONS 

For the discrete model of Eq. ([T]), in Sec. Ill Al we obtained exact moving solutions in the 
form of the sn Jacobi elliptic function. Solutions in the form of cn and dn Jacobi elliptic 
functions can also be constructed as well as the solutions having the form of 1/sn, 1/cn, 
and sndn/cn in analogy with [l^]. In Sec. Ill Bt from the sn solution, in the limit of m — > 1, 
we extracted the exact moving kink solution in the form of tanh, i.e., the corresponding 
hyperbolic function solution. 

Setting f = in the exact moving solutions obtained in this work one obtains the TI static 
solutions reported in 13]. In this sense, the results reported here generalize our previous 
results. We thus reveal the hidden connection between the static TI solutions and the exact 
moving solutions. Such solutions can be derived from a three-point map reducible to a 



two-point map (see Sec. Ill Fl) i.e., from an integrable map 37 1. 

We have demonstrated that the exact moving solutions to lattice and continuous equa- 
tions with competing nonlinear terms can have any large propagation velocity. Most of 
the high-speed solutions reported in the present study are solutions to the non-Hamiltonian 
variant of the considered 0^ models. However, the trigonometric solution described in Sec. 
Ill CI exists also in the Hamiltonian lattice and it also can have an arbitrary speed. 

While the problem of identifying traveling solutions in the one-dimensional context by 
now has a considerable literature associated with it, as evidenced above, identifying such 
solutions in higher dimensional problems is to a large extent an open question. While initial 
studies have demonstrated the possibility in some of these systems for traveling in both on- 



and off-lattice directions [38|, analytical results along the lines discussed here are essentially 
absent in that problem and could certainly assist in clarifying the potential of coherent 
structures for unhindered propagation in these higher dimensional settings. 
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